Search for: All records

We present a Bayesian framework for joint and coherent analyses of multimessenger binary neutron star signals. The method, implemented in our bajes infrastructure, incorporates a joint likelihood for multiple datasets, support for various semi-analytical kilonova models and numerical-relativity (NR) informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. As a first application, we analyze the gravitational-wave GW170817 and the kilonova AT2017gfo data. These results are then combined with the most recent X-ray pulsars analyses of PSR J0030+0451 and PSR J0740+6620 to obtain EOS constraints.Various constraints on the mass-radius diagram and neutron star properties are then obtained by resampling over a set of ten million parametrized EOS built under minimal assumptions. We find that a joint and coherent approach improves the inference of the extrinsic parameters (distance) and, among the instrinc parameters, the mass ratio. The inclusion of NR informed relations strongly improves over the case of using an agnostic prior on the intrinsic parameters. Comparing Bayes factors, we find that the two observations are better explained by the common source hypothesis only by assuming NR-informed relations. These relations break some of the degeneracies in the employed kN models. The EOS inference folding-in PSR J0952-0607 minimum-maximum mass, PSR J0030+0451 and PSR J0740+6620 data constrains, among other quantities, the neutron star radius to R1.4=12.30−0.56+0.81R1.4​=12.30−0.56+0.81​ km (R1.4=13.20−0.90+0.91R1.4​=13.20−0.90+0.91​ km) and the maximum mass to Mmax=2.28−0.17+0.25 M⊙Mmax​=2.28−0.17+0.25​ M⊙​ (Mmax=2.32−0.19+0.30 M⊙Mmax​=2.32−0.19+0.30​ M⊙​) where the ST+PDT (PDT-U) analysis of Vinciguerra et a (2023) for PSR J0030+0451 is employed. Hence, the systematics on PSR J0030+0451 data reduction currently dominate the mass-radius diagram constraints. 

Using a data set of approximately 2 million phenomenological equations of state consistent with observational constraints, we construct new equation-of-state-insensitive universal relations that exist between the multipolar tidal deformability parameters of neutron stars, Λl, for several high-order multipoles (l=5,6,7,8), and we consider finite-size effects of these high-order multipoles in waveform modeling. We also confirm the existence of a universal relation between the radius of the 1.4M⊙ NS, R1.4 and the reduced tidal parameter of the binary, Λ˜, and the chirp mass. We extend this relation to a large number of chirp masses and to the radii of isolated NSs of different mass M, RM. We find that there is an optimal value of M for every M such that the uncertainty in the estimate of RM is minimized when using the relation. We discuss the utility and implications of these relations for the upcoming LIGO O4 run and third-generation detectors.